Application of Pade interpolation to stationary state problems 
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If the small and large coupling behavior of a physical system can be computed perturbatively 
and expressed respectively as power series in a coupling parameter g and a Pade approximant 
embracing the two series can interpolate between these two limits and provide an accurate estimate 
of the system's behavior in the generally intractable intermediate coupling regime. The methodol- 
ogy and validity of this approach is illustrated by considering several stationary state problems in 
quantum mechanics. 



I. INTRODUCTION 

A Pade approximation is a formal transformation of 
the first n terms in the power series for a function f{x) 
into a rational function. The rational function R{x), 
known as the Pade approximant, typically comprises a 
ratio of two polynomials of a;, chosen in a way that a Tay- 
lor expansion of R{x) completely reproduces the original 
power series up to order n. In contrast to the truncated 
"mother" series, the Pade approximant is able to main- 
tain remarkable fidelity to the true f{x) for values of 
X well beyond the radius of convergence of the original 
power series, although how accurate it is, or how large 
I a; I may be before the approximation fails is impossible 
to establish in general. Nevertheless, its extraordinary 
predictive power has in the past been exploited in areas 
of quantum field theory and statistical physics. In quan- 
tum electrodynamics and quantum chromodynamics, for 
example, the Pade method has been shown to be an ef- 
fective means of both estimating unknown higher order 
terms as well as of summing the perturbation series for 
a physical observable that has been calculated to some 
finite order in the coupling constant.^ 

In a recent paper one of us has explored a different 
usage of the Pade approximation in which the Pade ap- 
proximant, constructed from the truncated power series 
of f{x) about two different points, serves to interpolate 
f{x) between the two points of expansion. In particu- 
lar, if one is able to compute perturbatively both the 
small X and large x behavior of f{x) (x > 0), and ex- 
press them respectively in a power series in x and a 
Pade approximant that simultaneously satisfies the two 
perturbation series will provide an accurate estimation 
of f{x) for the entire range of x, provided that f{x) is 
a sufficiently smooth function in this range. This ap- 
proach is known as the two-point Pade approximation,^ 
and is especially savoring from the perspective of duality 
in supersymmetric gauge theories)^ if one can infer per- 
turbatively the strong coupling behavior of a theory from 
its weakly coupled dual theory, a Pade interpolation of 
this limit with the perturbative weak coupling limit of the 
original theory will then give us a handle on the behavior 
of the theory for all coupling strengths. 

The Pade interpolation method may also be employed 
in some cases where there are no obvious expansion pa- 



rameters. An example is a system for which the Hamil- 
tonian H, although itself not exactly soluble, can be ex- 
pressed as a sum of two constituent Hamiltonians, Hi 
and H2, both with known exact solutions. It was pro- 
posed in Ref. 2 that, by tagging either Hi or H2 with an 
interpolation parameter A, the two required perturbative 
series in opposite limits of A could be generated, and 
an estimate for H would correspond to evaluating the 
Pade approximant with A — 1. This method was applied 
in Ref. 2 to heavy quarkonium systems with reasonable 
success. 

The quarkonium example demonstrates the power and 
usefulness of Pade interpolation for treating a certain 
class of stationary state problems. The technique may 
serve as an improvement or a supplement to perturbation 
theory typically taught in courses on quantum mechan- 
ics. The purpose of this note is to present a pedagogi- 
cal discourse of the methodology of Pade interpolation. 
For the examples considered below, the Pade interpola- 
tion method is shown to be stable to higher orders, and 
yield particularly good results when the usual perturba- 
tive method fails. 



II. THE METHOD 

Consider a system governed by a Hamiltonian H that 
has no known solution, but which may be expressed as a 
sum of two parts: H = Hi + H2, where Hi and H2 are 
individually exactly soluble. We are interested in finding 
the spectrum of H. In order to implement the Pade inter- 
polation for this purpose, we introduce the interpolation 
Hamiltonian 



H{X) = Hi + XH2 , 



(1) 



where the interpolation parameter A is real and positive. 
Suppose first that A ^ 1. We may then regard H2 as a 
perturbation, and calculate the eigenvalues of H as power 
series in A (the subscript j labels the eigenvalues): 



E<{X) = Co + ciA + C2A2 + . . . + c„,A™ . 
In the opposite A » 1 limit, we rewrite Eq. as 



iJ'(A) = A 



H2 



(2) 



(3) 
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and a second perturbative solution of H{X) follows from 
treating 1/A as a small parameter: 



E>iX) = X 



bo 



Mi 



+ b2 ( i 



(4) 

Note that both Hi and H2 must admit stationary states 
and be able to be treated as perturbations for the method 
to work. An example of a Hamiltonian that admits 
bound-state solutions but that cannot be treated as a 
perturbation is the Hamiltonian for an infinite rectangu- 
lar potential well. 

A generic Pade approximant for the energy eigenvalues 
has the form 



A- (A) 



N 



l + qiX + q2X'^ + 



qKlX'^ ' 



(5) 



where the N+AI +1 coefficients are determined by match- 
ing order by order the power series expansion of the Pade 
approximant with the perturbative results. For bg ^ 0, 
the nature of (A) in Eq. Q demands that the polyno- 
mials in the numerator and the denominator of the Pade 
approximant differ by one degree such that M — N — \. 
Furthermore, suppose that we solve ll(X) for small and 
large A to the same order in perturbation theory, that is, 
m = n. (This case is just an illustration and is not a 
necessary condition for implementing the Pade interpo- 
lation, although the accuracy of the approximation will 
depend on m and n.) Then Eqs. ij^Jl and (0J| together fur- 
nish 2n + 2 simultaneous equations for the 2N unknown 
coefficients p and q, and consequently the polynomials 
in the Pade approximant must satisfy N ~ n + 1 and 
M = n. The final step of setting A = 1 in the Pade ap- 
proximant yields an estimate for the eigenvalues of the 
original Hamiltonian H . 



III. EXAMPLES 

We demonstrate here the validity of the Pade interpo- 
lation method by way of two examples. Consider first a 
simple two-state system described by the Hamiltonian 



H 



ha. 



(6) 



where a and h are real parameters, and <Ty and are 
the Pauli matrices. For example, this Hamiltonian can 
represent the interaction energy of a charged spin-1/2 
particle in a magnetic field B = {0, By, Bz). In this case. 



ghBy 



and b - 



gfiB, 



(7) 



where g is the gyromagnetic ratio of the particle. This 
example is trivial in the sense that H can be easily diag- 
onalized to yield the exact eigenvalues 



E± = ±^Ja? + b'^ . 



(8) 



However, the comparison of these exact results with the 
approximate eigenvalues obtained below by Pade inter- 
polation will provide a way to gauge the accuracy of the 
approximation method. 

For \a\ ^ \b\, for example, if the magnetic field is al- 
most aligned with the y-axis, the baz term in H may be 
treated as a perturbation. To be more explicit, we may 
express the Hamiltonian as 



(9) 



which has the form of Eq. except that A is equal to 
b/a and corresponds to a physical expansion parameter. 
If we calculate the energy eigenvalues to second order in 
perturbation theory, we find 



E. 



a|>|fc| 



= ±|a| 1 



2a2 



(10) 



For \a\ ^ l&l, for example, if the magnetic field is almost 
parallel to the 2-axis, the aay term in H can be regarded 
as a perturbation. We find that, again to second order in 
perturbation theory (with a/b as the expansion parame- 
ter), the eigenvalues of H are now given by 



E. 



6|>|a| 



(11) 



A Pade approximant that interpolates these two limits 
of the energy eigenvalues can now be constructed. For 
the eigenvalue E^ , we find 



E^^ = 



|b/ap-t-f|Vap + ||b/a| + l 
\b/a\^ + l\b/a\ + l 



(12) 



This Pade approximant is uniquely determined from the 
perturbative expansions for E^ given in Eqs. H10() and 
(dH). Table d compares this Pade interpolation result 
with the exact eigenvalue, Eq. JHJ, for various values of 
the parameter \b/a\. We see that the Pade interpola- 
tion yields an approximation that is within 1% of the 
exact result for all values of \b/a\. This simple example 
demonstrates the potential power of the Pade interpola- 
tion technique: by simply computing the leading pertur- 
bative corrections for small and large \b/a\, one obtains 
a very accurate approximation to the eigenvalues for all 
values of \b/a\. 

As a second example, consider a single particle subject 
to a one-dimensional linear plus harmonic oscillator po- 
tential. (The Pade interpolation technique we shall use to 
solve this problem is similar to that applied in Ref. "5 to 
nonrelativistic quarkonium systems.) The Hamiltonian 
that describes the motion of the particle is 

where 7 > 0. (The latter condition is necessary in order 
for the Hamiltonian Hi in Eq. (|15|) to admit stationary 
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state solutions.) V{x) represents the rigid wall potential: 



V{x) = 







a; > 



oo . otherwise 



(14) 



The presence of V{x) restricts the particle's motion to 
be along the positive x-axis. 

The Hamiltonian in Eq. IjlSI) without V{x) appears in 
many textbooks on quantum mechanics,'' and can be eas- 
ily solved by transforming to a new coordinate x' with 
the origin at x = —^/(muP'). The presence of the rigid 
wall potential, however, requires all wave functions to 
vanish for a; < and renders such a coordinate redefini- 
tion useless. The Hamiltonian with 7 = is also a typical 
textbook problem^ that is exactly soluble; the boundary 
condition due to the rigid wall forces all energy eigenfunc- 
tions to vanish at the origin, which implies that only the 
harmonic oscillator states with odd parity are allowed. 

Solving the complete Hamiltonian H13|) is a somewhat 
more challenging task. In particular, if the linear and 
quadratic potentials are comparable, conventional per- 
turbative methods are not applicable. We shall therefore 
resort to the method outlined in Sec.mto find its eigen- 
values. Note that as long as we confine the particle's 
motion to the positive x branch, and impose the bound- 
ary condition that all eigenfunctions vanish at x = 0, we 
may drop the rigid wall potential V{x) in Eq. H13() . The 
resulting Hamiltonian can then be cast in the form of 
Eq. Q, with 



and 



Ho 



4 dx^ 

4 



■jx, 



r 



(15) 
(16) 



where, for simplicity, we have set h = 1, m = 1 and uj = 1 
such that the arbitrary parameter 7 alone regulates the 
relative importance of the two potential energy terms. 
Note that it is also necessary to split the kinetic energy 
term. Here, we have arbitrarily put half of the original 
kinetic energy term into each of the sub-Hamiltonians in 
Eqs. ((T3|l and As we will see (in the last paragraph 
of this section), better accuracy will generally be achieved 
in Pade interpolation if a larger fraction of the kinetic 
energy is included in the sub-Hamiltonian containing the 
dominant potential energy term. 

The solutions to Hi are the familiar Airy functions 
Ai{z), with z — {2/^)'^^^{-fx — e^); the energy eigenvalues 
are determined by the roots of Ai{z). On the other 
hand, H2 is solved by exp(— ^^/2)ft,j (^), where hj{^) are 
the Hermite polynomials of degree j, ^ = = 

(j -I- 1/2) /^/2 are the allowed energies, and the index j 
must be an odd integer in order to satisfy the boundary 
condition ■0^(0) = 0, where ipj{x) denotes the stationary 
state wave functions. 

We now proceed to perform the relevant perturbative 
calculations. We have evaluated to first, second, and 



third order in A and 1/A the approximate ground and 
first excited state energies for 7=1, that is, when the 
linear and quadratic potential energy terms are compa- 
rable, and have formed the unique Pade approximant for 
each instance. Because closed-form expressions for in- 
tegrals involving Airy functions generally do not exist, 
we did the exercise numerically. As an illustration, the 
analysis of the ground state generates the two series, 

E< = 1.47292 + 1.06950A- 0.0131354A2 + ..., (17) 

S^' = 1.06006A+ 1.47918 - 0.00467253^ + ... (18) 

A 

The first (second) order Pade approximant, 
1.47292 -f- 3.14779A + 1.49659A2 



E, 



(19) 



E, 



PA _ 

1 + 1.41100A 

PA 1.47292 + 5.36462A + 6.07765A2 + 2.14021A3 



0,2nd 



1 + 2.91607A + 2.01781A2 



(20) 



follows from manipulating the first two (three) terms of 
Eqs. (HZI and lO. Tables ITIl and lIlTl contain a summary 
of the results for the ground and first excited states. 

It is instructive to compare these results with those one 
would obtain from conventional perturbative calculations 
alone. Because there is no preference for either of the two 
potential energy terms, we consider both of the following 
parameterizations of the Hamiltonian: 



and 



where a and f3 are the small parameters that are even- 
tually set to 1. As seen in Tables IhI and Hill the Pade 
interpolation gives by far the most stable results. A fur- 
ther comparison with exact solutions from the numerical 
integration of the Schrodinger equation, also given in Ta- 
bleslllland lllll exemplifies the validity of the method. We 
have also checked the accuracy of the method for higher 
excited states. The approximate energies obtained, even 
to first order in the perturbation parameters, are always 
accurate to within 1% of their exact values. 

For completeness, we have examined situations in 
which one potential energy term is dominant, and per- 
turbative calculations on the smaller term alone are ex- 
pected to yield reasonably accurate results. This is cer- 
tainly the case. However, as seen in Tables IIVI and 
the Pade interpolation is able to do a better job, pro- 
vided that the original kinetic energy term is distributed 
among the two sub-Hamiltonians l|15|) and (|16|) in a way 
that reflects the relative significance of the two potential 
energy terms. We have also studied the effects of dis- 
tributing the kinetic energy unevenly between the two 
sub-Hamiltonians in the 7 = 1 case. As shown in Ta- 
bles I VII and IVIII rather good estimates of the exact re- 
sults can be achieved regardless of how the kinetic energy 



4 



is distributed, particularly if one goes to higher order. 
However, the best accuracy is obtained if somewhat less 
kinetic energy (40% to be precise) is included in Hi, es- 
pecially for the first excited state. This result can be un- 
derstood from Tables ITTl and IIIII which show that the per- 
turbation series for H'{X) converges faster than H{X) to 
the exact result. This behavior in turn suggests that for 
7 = 1, the linear potential is weaker than the quadratic 
potential. Hence, according to the results in Tables IIVI 
and a more accurate Pade approximant will be ob- 
tained by underweighting the kinetic energy in Hi. Un- 
fortunately, there are no quantitative rules for how the 
kinetic energy should be distributed among the two sub- 
Hamiltonians. Tables IvTl and IvTll suggest that a 50 : 50 
split should produce reasonably good estimates. 



Zeeman effects in hydrogen. Exact results for the n = 2 
level can be found in Ref. They involve square-root 
functions of the expansion parameter (the magnitude of 
the magnetic field) , similar to the first example discussed 
in Sec. IIIII See also Ref. which discusses the case of the 
two-dimensional hydrogen atom. 

The use of Pade interpolation is of course not limited to 
quantum mechanical problems, because all that is needed 
is an expansion parameter, be it a physical one as in 
Eq. or an artificial one such as the interpolation pa- 
rameter in Eq. Q , for which the behavior of the physical 
system can be calculated or measured when the parame- 
ter is small as well as when it is large. We encourage the 
reader to find other applications of this useful approxi- 
mation scheme. 



IV. CONCLUSION 

The stationary state problems considered here provide 
a good illustration of the power of Pade interpolation for 
problems for which exact solutions are difficult to obtain 
and ordinary perturbation methods are not applicable. 
For practice, the interested reader may wish to apply 
the method to interpolate the strong-field and weak-field 
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TABLE I: The eigenvalue E+ (in units of |a|) of the two-state 
system for various values of the parameter \b/a\. 



b/a\ Pade Exact 

0.1 1.00517 1.00499 

0.5 1.12500 1.11803 

1 1.42857 1.41421 

2 2.25000 2.23607 
10 10.0517 10.0499 
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TABLE II: Ground state energy, 7 = 1. Columns two to 
five display the first, second, and third order perturbative 
solutions to the Hamiltonians H{a), H{l3), H{X), and H'{X), 
where the perturbation parameters a, {3, A, and 1/A are all 
set to unity. Results from the Pade interpolation of H{\) 
appear in column six. These are to be compared with the 
exact energy, shown in the bottom, obtained from numerical 
integration of the Schrodingor equation. 

H{a) HjP) H{X) H'jX) Pade 

1st order 2.62838 2.77411 2.54242 2.53984 2.53724 

2nd order 2.51908 2.30374 2.52928 2.53517 2.53720 

3rd order 2.54121 2.88137 2.54998 2.53882 2.53720 



TABLE III: First excited state energy, 7 = 1. See Table |n| 
caption for a detailed description. 

H{a) H{(3) H{\) H'{\) Pade 

1st order 5.19257 6.05194 5.20217 5.13559 5.10483 

2nd order 5.09417 3.67756 4.90789 5.08365 5.10380 

3rd order 5.08881 8.19137 5.53588 5.11655 5.10333 

Exact 5.10382 



Exact 



2.53720 
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TABLE IV: Perturbation theory versus Pade interpolation: 
dominant quadratic potential, 7 = 0.1. Columns two and 
three show respectively the results from standard perturbative 
calculations with the linear potential as the perturbation, and 
the corresponding estimates from the Pade interpolation of 
H{X), where the original kinetic energy term is distributed 
among the sub-Hamiltonians 1151 and in the ratio 1 : 9. 
The exact energies, obtained from numerically integrating the 
Schrodinger equation, are also displayed. 

Perturbation Theory Pade 

Ground state 

1st order 1.61284 1.61178 

2nd order 1.61174 1.61177 

3rd order 1.61177 1.61177 



TABLE V: Perturbation theory versus Pade interpolation: 
dominant linear potential, 7 — 10. Columns two and three 
show respectively the results from standard perturbative cal- 
culations with the quadratic potential as the perturbation, 
and the corresponding estimates from Pade interpolation of 
H{X), where the original kinetic energy term is distributed 
among the sub-Hamiltonians l|15^ and l|16|l in the ratio 9:1. 
The exact energies, obtained from numerically integrating the 
Schrodinger equation, are also displayed. 

Perturbation Theory Pade 

Ground state 

1st order 8.81152 8.80704 

2nd order 8.80681 8.80706 

3rd order 8.80708 8.80706 



Exact 



1.61177 



First excited state 
1st order 
2nd order 
3rd order 



3.66926 
3.66827 
3.66827 



3.66844 
3.66828 
3.66828 



Exact 



8.80706 



First excited state 
1st order 
2nd order 
3rd order 



15.6650 
15.6412 
15.6433 



15.6432 
15.6431 
15.6431 



Exact 



3.66828 



Exact 



15.6431 
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TABLE VI: Comparisons similar to the last three columns 
in Table UTI except the original kinetic energy is distributed 
among the sub-Hamiltonians 1151 and 1161 1 in the ratios indi- 
cated. 



1 : 9 

1st order 
2nd order 
3rd order 



2 : 8 

1st order 
2nd order 
3rd order 



3 : 7 

1st order 
2nd order 
3rd order 



4 : 6 

1st order 
2nd order 
3rd order 



5 : 5 

1st order 
2nd order 
3rd order 



1st order 
2nd order 
3rd order 



7 : 3 

1st order 
2nd order 
3rd order 



H{X) 



3.64332 
-2.96766 
33.5854 



2.84633 
1.66075 
5.07816 



2.62009 
2.35777 
2.87926 



2.54956 
2.51819 
2.56636 



2.54242 
2.52928 
2.54998 



2.56633 
2.49717 
2.59773 



2.60713 
2.45110 
2.66226 



H'iX) 



2.60113 
2.52016 
2.54173 



2.57650 
2.52341 
2.54187 



2.55603 
2.52850 
2.54110 



2.54229 
2.53396 
2.53924 



2.53984 
2.53517 
2.53882 



2.55756 
2.51796 
2.55661 



2.61519 
2.43389 
2.68987 



Fade 



2.54126 
2.53730 
2.53721 



2.53717 
2.53736 
2.53720 



2.53703 
2.53720 
2.53720 



2.53718 
2.53720 
2.53720 



2.53724 
2.53720 
2.53720 



2.52772 
2.53718 
2.53720 



2.54018 
2.53704 
2.53721 



TABLE VII: Comparisons similar to the last three columns 
in Table UTTl except the original kinetic energy is distributed 
among the sub-Hamiltonians 1151 and 1161 1 in the ratios indi- 
cated. 



1 : 9 

1st order 
2nd order 
3rd order 



2 : 8 

1st order 
2nd order 
3rd order 



3 : 7 

1st order 
2nd order 
3rd order 



4 : 6 

1st order 
2nd order 
3rd order 



5 : 5 

1st order 
2nd order 
3rd order 



6 : 4 

1st order 
2nd order 
3rd order 



7 : 3 

1st order 
2nd order 
3rd order 



HiX) 



6.62887 
-3.53137 
80.5535 



5.38749 
4.10020 
10.1569 



5.11949 
5.05046 
5.40290 



5.11002 
5.08677 
5.14789 



5.20217 
4.90789 
5.53588 



5.34182 
4.66966 
6.03944 



5.50568 
4.41718 
6.56543 



H'iX) 



5.15343 
5.09570 
5.09344 



5.12254 
5.09973 
5.09757 



5.10398 
5.10423 
5.10063 



5.10443 
5.10370 
5.10373 



5.13559 
5.08365 
5.11655 



5.21984 
5.00364 
5.18832 



5.40620 
4.73387 
5.57852 



Fade 



5.10321 
5.10468 
5.10330 



5.10284 
5.10380 
5.10302 



5.10341 
5.10382 
5.10344 



5.10377 
5.10382 
5.10382 



5.10483 
5.10380 
5.10333 



5.10989 
5.10355 
5.04811 



5.12506 
5.10225 
5.10083 



1st order 
2nd order 
3rd order 



9 : 1 

1st order 
2nd order 
3rd order 



2.65771 
2.40135 
2.73243 



2.71414 
2.35175 
2.80562 



2.76705 
2.07070 
3.65519 



3.24341 
-0.240897 
16.5922 



2.54828 
2.53623 
2.53729 



2.57328 
2.53211 
2.53797 



1st order 
2nd order 
3rd order 



9 : 1 

1st order 
2nd order 
3rd order 



5.68231 
4.16511 
7.09878 



5.86556 
3.91817 
7.63987 



5.82763 
3.73063 
8.09921 



7.03919 
-2.12398 
39.8902 



5.16388 
5.09690 
5.10182 



5.26810 
5.07466 
5.10548 



